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The interplay of strong and weak decay amplitudes for B — > ■kttK and B — > KKK, with 
the TTTT and KK pairs interacting in isospin-0 S-wave, is analyzed for tttt effective mass 
from threshold to 1.2 GeV. To improve agreement with experiment of a factorization 
approach with some QCD corrections, addition of long-distance contributions, called 
charming penguins is necessary. 



1. Introduction 

Hadronic two- and three-body B-meson weak decays are a rich source of informa- 
tion to study CP violation within the Standard Model (and beyond). On Dalitz 
plots, many resonances are distinguishable but with entangled interference patterns. 
Moreover the understanding of the final state interactions is weighty to attain accu- 
rate values of the Cabibbo-Kobayashi-Maskawa (CKM) matrix elements. We report 
here on our recent studj^ (FKLL) of i? ^ ttttK and B KKK decays. We con- 
sider the TTTT or KK pairs, denoted by (7r7r)5 or {KK)s, to interact in isospin-0 
5- wave in the mTTTr effective mass range from threshold to 1.2 GeV. We extend to 
the /o(980) case the approach of Gardner and Meil3neiEl (GM) who have examined 
the effect of the /o(600) (or a) on the i?° 7r+7r~7r° decays. The decay ampli- 
tudes are calculated within the QGD factorization approach. Due to cancellations 
between penguin amplitudes, the B /o(980)i4r branching fractions are too small 
compared to data. We add long-distance contributions, called charming penguins. 
Expressions for the B decay amplitudes and description of the final state inter- 
actions, are given in section II. Comparison of our results to the data and some 
conclusions are presented in section III. 
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2. B decay amplitudes and final state interactions 

From the effective weak Hamiltonian, H^ff, and the operator product expansion, 
the B decay anipHtudc into mesons Ml and M2 iP (/ stands for quarks d or s) 



q—u.c 
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(1) 



In (1) Gp = 1.66 X 10"^ GeV-2 is the Fermi constant, V^/^^ represents the CKM 
matrix-element factors and n is a. renormalization scale. The perturbative Wilson 
coefficients Ck{n) characterize the short-distance physics. The non-perturbative 



hadronic matrix elements (AfiM2|0'''(/i)|i?) describe the long-distance physics. 



The local operators Ol'^fi) "effectively" govern the decay. The main task of the 
theory is to compute these matrix elements in a reliable way. Beneke, Buchalla, 
Neubert and SachrajdsPhave developed the following QCD factorization formula, 



{MiM2\0l\^i)\B) ^ (A/i|Ji|0)(M2|J2|B) 



1 + ^ r„< -f 0{KQCDlmi,) 



(2) 



In (2) Ji and J2 are bilinear quark currents, (Mi|Ji|0) is the meson Mi decay 
constant and {M2\J2\B) represents the transition form factor of the meson B into 
the meson A/2. The radiative corrections r„ have been calculated to order n = 1, 
tts being the strong coupling constant. The term 0{KQCD/'mh) in Eq. (0) repre- 
sents power corrections and if set to zero, together with r„, one recovers the naive 
factorization formula as appHed by Bauer, Stech and WirbeP. 

The possible quark line diagrams for the B^ — > (7r7r)5A'~ and B^ {KK)sK^ 
are shown in Fig. 1 of FKL1J2. For the 5° decays there is no tree diagram (a) , only 
penguin diagrams similar to (b) or (c) diagrams exist. Within the QCD factorization 
framework, just reminded above, the B~ {tt^tt~)sK~ decay amplitude is 

{{7:+7:-)sK-\H,ff\B-) = ^^i^x[P{m^n){Ua + Ut) + C{m^^)]Trim^.) + 

)}• (3) 

The constant x, to be fitted to reproduce the averaged branching fraction B\B^ 
/o(980)is:=^], is estimated to 30 GeV^^ from the /o(980) decay propertie!*li. The 
non-strange and strange pion scalar form factors r"(TO7r7r) and rf(m7r7r) describe 
the final-state interactions. In (3) P^rriT^T^) corresponds to the product of the kaon 
decay constant fK times the transition form factor of B into (7r7r)5, F^^^''^'^ (M^), 



P{m^^) = f^{M'^ - mi^)Fr'^->^{Mf,). (4) 

The function Q(mTr7r) is proportional to the product of the (7r7r)s " de cay c onstant" 
times the transition form factor of B into K, F^^^ [rn^^), and readiil^, 

g(™,,) = (Ml - Ml) Fr^" {ml.) , (5) 

vib — nis 



3 



with Bo = IP. = <l(jnu + md). 

The functions Ua^ Ub and Uc, given in terms of the coefficient a I I I I at 
of the decay operators, are calculated from the Wilson coefficients. They cor- 
respond to the contributions of Figs. 1(a), 1(b) and 1(c) of FKLlff, respec- 
tively. One has Ua = VubV:,ai, Ub = KbK*. W - a^) + ^cfcK*s (a^ " agr) = 
VubV:, K + (ag -<)r]+ VtbV,* (agr - a^) and = -KbK*.^ " VcbV^Xe = 
VubV*g (ag - flg) -I- VtbVtlag. We have used the unitarity of the CKM matrix el- 
ements, VubV*g + VcbV*g + VtbVf*g = 0. In the reduction of the four-quark opera- 
tors, split into product of two- matrix elements (Eq. (0)), a Fierz transformation 
is needed for the penguin diagrams in order to match the flavor quantum num- 
bers of the quark currents to those of the hadronJSl The chiral factor r is equal to 
2M'^/[(mb+mu){ms+mu)]. Here and in Eq. jSJ, quark masses appear as one makes 
use of the Dirac equation. The ai and o^'q coefficients at next-to-leading order in as 
are given m Eq. (35) of Beneke and NeuberiP (BN). In the numerical calculations 
we take the values of table III of de Groot, Cottingham and Whittinghan:^^ (GCW) 

at the scale u = 2.1 GeV. We do not include the small contributions from hard 
. fTI . . . 

spectator interactions , annihilation and clcctroweak diagrams. 

To improve their fit to charmless hadronic two-body i? decays with pseudoscalar- 
pseudoscalar and pseudoscalar-vector decays, GCW^Sl, following Ciuchini et al.l^, 
introduce the charming penguin term C(rn), the parametrization of which is re- 
minded in Eq. (6) of FKLlW Expressions for the other B {tt'k)sK and for the 
B {KK)sK amplitudes are given in FKLlP. If one neglects the penguin con- 
traction terms (see Eq. (35), (39) and (41) of BjJ^) then a^Q — al^ — a^ Q and our 
amplitude Q is similar to that of the B^ cttt" deca}El(see their Eq. (25)). Note 
the near cancellatiorP in the b —* s transition of the two penguin contributions as 
r ~ 1 and 04 ~ Og. 

The r"'''(m) obey unitary tttt and KK coupled channel equationiM We use the 
solution A of the two-body scattering matrix of Kamihski, Lesniak and Loisear^l 
Below the KK threshold r"'^*(m) = i?"''(m) cos (m)e*'^'"'("). |r"'"| is maxi- 
mum for (Jtttt — 180°, which corresponds to the /o(980) production. The i?"'''(m), 
depicting the for mat ion of mesons prior to rescattering, have been obtained by 
Meifiner and OUeiEl. 



3. Results and conclusions 

Our results with the charming-penguin C(m) terms of GCV^, model I, and of 
Ciuchini, Franco, MartineUi, Masicro, Picrini and SilvestrinP^ (CFMMPS), model 
II, are compared to the dat^^'^-i—--. i^ Table 1. Values of all parameters used 
are given m FKLlJU. The effective m,r7r mass distribution for B^ tt+tt K^ of 
our model II is compared with Belle^ in Fig. 1. The sharp maximum near 1 
GeV is from the /o(980). Near 0.5 GeV the broad bump is related to the /o(600). 
Events near 0.8 GeV are due to the B^ p^K^ decay not considered in our 
model. Similar agreement with m^rTr effective mass distributions of BaBar and Belle 
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Figure 1. Comparison of our model II with Belle dati 



Table 1. Average branching fractions B in units of 10 ^, direct asymmetries A.cp, -4 and 
S parameters of our model compared to d ata. Model errors come from C(m) parameters 
fitted by GCV\^ model I, and CFMMP^ model II. For a qualitative comparison we give in 
parenthesis the experimental B — > KKK decay observables obtained by considerin g the full 
range of the KK invariant mass rather than the upper limilB of 1.1 GeV. For dat J^'^ l l ^ we 
only quote statistical errors. 



B decay mode 




Average 
HFAG's valuesis 


Model I 
X = 33.5 GeV-i 


Model II 
X = 23.5 GeV-i 


B± ^/o{980)K±, 


B 




8.49 (fit) 


8.46 (fit) 


/O — > TT^TT^ 


Acp 


-0.02 ± 0.07 


-0.52 ± 0.12 


0.20 ±0.20 


^ /o(980)i<'0, 


B 


6.0 ± 1.6 


5.9 ± 1.6 


5.8 ± 2.8 




A 


-0.14 ±0.22 


0.01 ±0.10 


0.0004 ± 0.0010 




S 


-0.39 ± 0.26 


-0.63 ± 0.09 


-0.77 ± 0.0004 


B± ^ {K+K-)gK^ 


B 

Acp 


< 2.913 

no data 


1.8 ±0.4 
-0.44 ±0.12 


1.7±0.7 
0.29 ±0.21 




B 


no data 


1.1 ±0.3 


1.2 ±0.5 


B^^{K+K-)^K% 


A 

S 


(-0.09 ±0.10) 
(-0.55 ±0.22)14 
(-0.74 ±0.27)18 


0.01 ±0.10 
-0.64 ±0.09 


0.001 ±0.001 
-0.77 ±0.0006 




B 


no data 


1.1 ±0.3 


1.2 ±0.5 




A 


(0.41 ±0.21) 


0.01 ±0.10 


0.001 ±0.001 




S 


(-0.26 ±0.34) 


-0.64 ±0.09 


-0.77 ±0.0006 



for charged and neutral B decays can be seen in the figures shown in 

To conclude, our effective mass distributions and branching fractions compare 
well with data. The long distance contribution C(m) is important. If C{m) = 0, 
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B{B° TT+n-K'^) is too small by a factor - 18 and B{B^ ■n+TT-R^) by a 
factor ~ 4. If charming penguins are included the agreement with data is go od 
with a X compatible with the /o(980) properties. The recent determinatioii^ of 
Acp = -0.02 ± 0.07 for fo{980)K^, fo n+ir-, should help to determine 

better the charming-penguin C{m) parameters. In progress we are adding the P- 
wave contribution which will allow us to describe the p-meson region. 
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